Can TeVeS avoid Dark Matter on galactic scales? 
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A fully relativistic analysis of gravitational lensing in TeVeS is presented. By estimating the 
lensing masses for a set of six lenses from the CASTLES database, and then comparing them to 
the stellar mass, the deficit between the two is obtained and analysed. Considering a parametrised 
range for the TeVeS function fi(y), which controls the strength of the modification to gravity, it is 
found that on galactic scales TeVeS requires additional dark matter with the commonly used (J,(y). 
A soft dependence of the results on the cosmological framework and the TeVeS free parameters is 
discussed. For one particular form of fi(y), TeVeS is found to require very little dark matter. This 
choice is however ruled out by rotation curve data. The inability to simultaneously fit lensing and 
rotation curves for a single form of fi(y) is a challenge to a no dark matter TeVeS proposal. 

PACS numbers: 



The standard (ACDM) cosmological paradigm is based 
upon Cold Dark Matter (CDM), a cosmological constant 
A and classical general relativity/Friedman-Robertson- 
Lemaitre- Walker cosmology. Despite its enormous suc- 
cess and consistency with a plethora of astrophysical 
data, competing models have been proposed for the pri- 
mary reason of the still unknown nature of the dark en- 
ergy component and the current undetectability of dark 
matter. To explain the observed flat rotation curves 
of galaxies without Dark Matter, Milgrom [l[ proposed 
MOdified Newtonian Dynamics (MOND), based upon 
the relation f(\a\/ao)S = — V$n between the accelera- 
tion a and the Newtonian gravitational field <&n- The 
constant ao w 1.2 x 10~ 10 m/s 2 is motivated by the accel- 
eration found in the outer regions of galaxies where the 
rotation curve is flat. When /, assumed to be a positive, 
smooth, monotonic function, equals unity, usual Newto- 
nian dynamics hold, while when it approximately equals 
its argument, the deep MONDian regime applies. 

MOND has been successful in explaining the dynamics 
of disk galaxies, however it is less successful for clusters 
of galaxies. It was promoted [2j to a classical relativistic 
field theory by introducing a TEnsor, VEctor and Scalar 
field (TeVeS). TeVeS has been criticised as lacking a fun- 
damental theoretical motivation. Recently, it has been 
argued 0] that such a theory can emerge naturally within 
some string theory models. 

In Ref. |4j, where the lensing mass in MOND was 
compared to the stellar mass content of the lenses, it 
was found that comparable amounts of dark matter were 
needed in MOND to that required in the standard lens- 
ing scenario. This result is in contrast with attempts to 
explain the lensing data on galactic-cluster scales by in- 
troducing a 2 eV neutrino. In fact, this component of 
dark matter has been shown [f| to cluster on Mpc scales 
but not on galactic scales where the previous analysis was 
conducted. It was concluded that either lensing must op- 
erate in a qualitatively different way within the covariant 
"parent" theory of MOND, such as the TeVeS models, or 
dark matter should be considered within MOND even on 



galactic scales. 

In this letter TeVeS is examined in a similar way, 
namely by deriving the modified lensing equation and 
solving it numerically, to investigate whether it also 
shows a dark component from lensing. The (weak) de- 
pendence of the results on the cosmological models and 
TeVeS free parameters is then discussed. TeVeS [2.j is a 
bi-metric model in which matter and radiation does not 
feel the Einstein metric, g a p, appearing in the canon- 
ical kinetic term in the (effective) action, but a modi- 
fied "physical" metric, g a p, related to the Einstein metric 
by g af 3 = e-^g af3 - U a U p {e 2 ^ - e~ 2 *), where U„,<t> de- 
note the TeVeS vector and scalar field, respectively. The 
TeVeS action is: 
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where k, K are the coupling constants for the scalar, 
vector field, respectively; £ is a free scale length re- 
lated to a (c.f below, after Eq.©); a is an addi- 
tional non-dynamical scalar field; = U^ v — U v ^; 
A is a Lagrange multiplier implementing the constraint 
g a ^U a Uf3 — —1, which is completely fixed by variation 
of the action; the function F(kGa) is chosen to give 
the correct non-relativistic MONDian limit, with G re- 
lated to the Newtonian gravitational constant, Gn, by 
G = Gjv/(1 + -^j^ — 20 c ), where 4> c is the present day 
cosmological value of the scalar field. Covariant deriva- 
tives denoted by are taken with respect to g^ and 
indices are raised/lowered with the metric g^. A new 
function fi(y) is introduced as Q 
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Possible choices for the fi function will be discussed later. 

The isotropic spherically symmetric Einstein metric 
can be generically written as g al 3dx a dx 13 — —e v dt 2 + 



(dr 2 
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: ) ; both v and £ are functions 



of r. The physical metric g a p has the same form, with e v 
and e^, related to the Einstein metric functions through 
v = v + 2(j), C = C — 20. Isotropy implies <f> = 4>(r). As- 
suming an ideal pressureless matter fluid, T a p — pu a Uf3. 
Motivated by a homogeneous and isotropic cosmology, 
the vector field is considered time-like. The normalisation 
condition imposed by the Lagrange multiplier in Eq. ([1]) 
gives U a = {e- v ' 2 , 0,0,0). 

Previous attempts lfl, U , EH at lensing analysis in 
TeVeS have remained non-relativistic, considering only 
the effect of adding a scalar potential to the standard 
Newtonian potential. Hence they are insensitive to any 
unique features of TeVeS as a fully relativistic field the- 
ory. The analysis given here is used to solve (numerically) 
the TeVeS equations of motion and obtain explicit ex- 
pressions for the physical metric quantities i>, £ and thus 
derive the modified Birkhoff 's theorem for the TeVeS the- 
ory. These functions are used to obtain the deflection an- 
gle and find the lensing mass, which is then compared to 
the stellar mass content. Assuming a mass density profile 
m s (< r), within a radial distance r, leads to a system of 
differential equations to determine C and v. A transfor- 
mation to v, (, then gives the physical metric. The (tt) 
and (69) differential equations are: 
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where the prime denotes derivatives with re- 
spect to r and the mass profile is m s (< r) = 
47r / r pexp f | + ^ - 2</>r 2 ^ dr. The deflection an- 
gle reads @: 
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of the galaxy contained within the virial radius, r v ; r ; it 
also specifies the density profile. 

Making the approximation to s (< r) « M(< r), 
shown [3] to be correct to leading order of r, Eq. © is 
numerically solved and through a transformation to the 
physical metric, the deflection angle Eq. Q is obtained. 

A specification for the function a (see, Eq. [2]) is re- 
quired. Previous lensing studies [l2j on the non- 
relativistic scalar potential approach to TeVeS, adopted 
a model for n(y) given in Ref. 0]. It was however 
noted 14 , lal that when this choice is converted into its 



MONDian equivalent, rotation curve data are not well 
fitted. It was hence proposed to consider [l6} the MOND 
function which best fits the rotation curve data and then 
convert it into its TeVeS fJ-(y) analogue. Any intermedi- 
ate choice for p(y) can be parametrised by a. In what 
follows the results from testing this parametrised fJ-(y) in 
the fully relativistic treatment of lensing in TeVeS are 
presented in order to see if TeVeS can both fit rotation 
curve data and lensing data for a single choice of /J,(y). 
For TeVeS fi(y), F(ji) and MOND f(x) we use: 
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2 C(r )-p(r ) r 2 jg the observable impact parameter ployed m Refs 



where the parameter range is < a < 1 [16j . The a = 
case for fi(y) gives the weak and intermediate gravity 
limit of n(y) and F(/j,) is taken from its explicit form 0]. 
The a = 1 gives the function which better fits rotation 
curve data. Since the functions increase monotonically, 
the analysis could be confined to the extremes of the 
parameter space, i.e. the a — and a = 1 cases. The 
TeVeS parameters are [Io| k = 0.01, K = 0.01, I = 

\/fc6/(47rSao), (j) c = 0.001, 5=1 + K/2 - 20 c . 

The cosmological model used is (fi m! ^A,^k) = 
(0.3, 0.7, 0), though the effect of other choices is also con- 
sidered. The parameter b is fixed from the limit of the 
y(n) function as fx becomes << 1, w 6/z 2 . For our 
parametrised choice of n(y) b — 3, specifying then I. The 
deflection angle for a model system is then calculated 
with this choice of parameters. The results for the deflec- 
tion angle in TeVeS are compared to the deflection angle 
resulting from following the scalar potential method em- 

n Elm, as wen as gr and mond. 



and ro is the point of closest approach for the light ray. 

The Navarro-Frenk-White (NFW) profile jQ 
is used, which in Schwarzschild radial coordi- 
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The deflection angle results are shown in Fig. [T] 

To calculate a possible dark matter component, a sam- 
ple of double lensing systems from the Castles database 
is analysed and the mass of the lensing galaxy in GR, 
MOND, and TeVeS is calculated. By comparing the mass 
from lensing to the stellar mass content calculated from a 
comparison of photometry and stellar population synthc- 
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FIG. 1: Deflection angle curves for a NFW profile system 
(C = 10). The parameters of the system are Di = 1000 Mpc 
(the luminosity distance of the lens), M = 10 1 Mq, k — 
0.01, K = 0.01, <j> c = 0.001 and r vir = 15 Kpc. Left panel: 
a — case. Right panel: a = 1 case. The deflection angles 
decrease in MOND and TeVeS as a increases. 



lar masses (excluding one outlier), it is found that even 
in TeVeS on average the dark matter content is 48.5% 
when a = 0, and 34.3% when a — 1. The outlier lensing 
system BRI0952-0115 shows a mass overshoot of « 80% 
in TeVeS, i.e. the lensing mass in TeVeS is less than the 
stellar mass content. It is possible that this system is 
affected by some unknown lens environment effects such 
as an unseen cluster mass contribution as has been sug- 
gested of other lenses [111], though there is no data to 
conclude this at present. Overall the analysis shows that 
TeVeS finds it hard to explain the lensing observed in 
these systems using only the stellar content of the galax- 
ies, a conclusion which stands in contrast to that given 
in Refs. Ill, ll2| despite the two methods predicting sim- 
ilar deflection angles (see, FigO}. This is a problem for 
TeVeS, which will be made more explicit later on in this 
letter, when the analysis of varying the free parameters 
will be performed. 



sis using a Chabrier Initial Mass Function (IMF), as in 
Ref. [l3|, the mass deficit which belonged to the "dark" 
sector is found. Note that while the authors of Ref. @ 
claim that the stellar mass estimates are sensitive to the 
IMF used, it has been argued in Ref. Q that the main 
competitor of the Chabrier IMF, the Salpeter IMF, ap- 
pears to fit the data worse. Other realistic IMF choices 
differ from the Chabrier by an insignificant factor 
These arguments support the validity of the method em- 
ployed here. Finally, inverse ray tracing is used to calcu- 
late the mass Q • The mass estimates we obtain this way 
are shown in Table 1. 
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TABLE I: Mass estimates (in 1O 1O M units) for ACDM cos- 
mology: (f7 m ,fi A ,^k) = (0.3,0.7,0), k — K — 0.01, and a 
NFW halo profile with C = 10. Two different cases (a = 0, 1) 
for fi(y) parametrisation are considered. 

Figure [5] compares the mass estimates between GR, 
MOND and TeVeS for the two cases of n(y), with the 
mass difference being given as a function of the mass cal- 
culated using GR (left panels) and Ri ens /R e (right pan- 
els), where i?i on s is the distance out to which our mass 
estimates and the stellar mass estimates are calculated 
and R c is the half light radius, both given in Ref. [HI]. 
The top two panels are for a = 0, the bottom two are for 
a = 1. Comparing the lensing masses against the stel- 
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FIG. 2: A comparison of the need for dark matter for the 
lenses in GR, MOND and TeVeS, for a = (top panels) 
and a = 1 (bottom panels). The right hand panels give the 
dark matter requirements as a function of the ratio Ri en s/Re- 
Any dark matter needed in TeVeS is especially significant for 
those masses which probe further out, where the modification 
to gravity is larger. 

The results show that the case a = 1, which is a choice 
specifically adopted to fit the rotation curve data, re- 
quires an even larger amount of dark matter. Thus, at- 
tempts to fit TeVeS to a no dark matter scenario using 
the freedom in the parameters of the class of /1 functions 
used here and in the literature so far, would imply that 
the theory fits poorly the rotation curves data, which was 
the original motivation for modifying the gravitational 
behaviour. 

The effect of different cosmological parameters has also 
been examined to see how they alter the results for the 
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a = case, corresponding to the minimum amount of 
required dark matter. It is found that the lensing re- 
sults are insensitive to the precise cosmological param- 
eters, which is not surprising, since the observational 
constraints mostly impose limits on the luminosity and 
angular distance scales For completeness we state 
the results, all of which point towards fluctuations in 
the amount of dark matter well within the error lim- 
its. In particular, for the case [T3| (Cl m , Qk) = 
(0.03, 0.46, 0.51), averaging over the six galaxies, we find 
a 5.7% increase in the amount of dark matter required. 
For the case (f2 m , f^A, ^k) = (0.23, 0.78, 0), considered in 
Ref. to fit the Cosmic Microwave Background data 
with TeVeS, one finds a corresponding average decrease 
by 2.2 %, while for <f> c = 0.01, one finds an average de- 
crease by 1.1% . 

Finally, the fc, K parameters are varied independently 
to check on the robustness of our claims. In particular, we 
examine the lensing system HE1104-1805, which within 
the TeVeS approach (with a = 0) requires the largest 
amount of dark matter. We find that the variation of the 
parameter k has considerably smaller effects than that of 
K. This supports the dominant role played by the vector 
field in TeVeS. A similar result has also been pointed out 
in Ref. [17| . but from a different perspective. Within the 
allowed parameter space 1 < K < 10~ 5 , 1 < k < 10~ 5 , 
implied by rotation curve data and solar system tests 
of gravity, the amount of dark matter required in the 
system HE1104-1805 is negligible only for values of K 
higher than 0.1, which however is excluded by gravita- 
tional measurements at solar-system scales. 
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FIG. 3: The dark matter content of the lenses for a = — 1 in 
the parametrisation of f(x) in MOND and fi(y) in TeVeS. 

We also considered other values of a, outside the region 
a G [0, 1], considered in the literature so far. We found 
that the a = — 1 case appears not to require substan- 
tial amount of dark matter to explain the gravitational 
lensing, and hence it could provide an example of an al- 
ternative to the dark matter scenario. However, such 
a model is in conflict with data from rotation curves of 



galaxies. 

In conclusion, the above results show a soft dependence 
on the free parameters of TeVeS and the cosmological 
model adopted, but a rather strong one on the form of 
the [i(y) function. Our analysis in this letter points to- 
wards the fact that TeVeS, at least within the class of 
models considered so far in the literature, cannot sur- 
vive both gravitational lensing and rotation curve tests. 
However, we cannot yet exclude completely the possibil- 
ity, admittedly remote, that a class of fi functions, or 
more complicated equivalents thereof, can be found such 
that alternative to dark matter scenarios are at play. 
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